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OBJECTIVES 

(a)  The  long  range  goal  of  the  research  is  to  provide  an  accurate  parameteri- 
zation of  the  acoustic  field  scattered  by  a randomly  rough  surface.  The 
research  is  designed  to  evaluate  existing  rough  surface  scattering  models  in 
view  of  their  application  to  scattering  from  the  under  surface  of  the  ice. 

(b)  The  immediate  objectives  of  this  problem  are  to  construct  random  surfaces 
with  surface  height  distributions  that  are  measurably  Gaussian  and  have  known 
correlation  lengths;  display  the  acoustic  field  scattered  from  the  surface 
by  means  of  schlieren  visualization;  measure  the  scattered  field  with  under- 
water transducers;  and  compare  the  results  to  existing  models  for  scattering 
of  sound  by  random  surfaces. 

ACTIVITIES/ACCOMPLISHMENTS 

The  investigations  undertaken  through  both  the  task  at  the  Naval  Research 
Laboratory  and  the  contract  at  The  American  University  have  been  successfully 
completed. 

Two  rough  surface  topographies  were  constructed  under  computer  control 
to  have  surface  height  distributions  and  slope  distributions  that  differed 
from  Gaussian  statistics  by  no  more  than  0.0036  (mean  squared  error).  Table  I 
lists  the  principal  surface  parameters  for  both  surfaces.  Rough  surface  I is 
shown  in  Fig.  1. 

Each  surface  was  used  in  a series  of  underwater  acoustic  scattering 
measurements  which  employed  three  different  experimental  techniques  to  deter- 
mine the  scattering  properties  of  the  surface.  The  geometry  diagrammed  in 
Fig.  2 was  maintained  throughout  all  experiments. 


Schlieren  Visualization. 


A typical  schlieren  visualization  system  is  shown  in  Fig.  3.  This 
method  was  used  to  provide  a qualitative  view  of  the  entire  scattered  field. 

The  photographs  in  Fig.  4 show  a 1 MHz  pulse  directed  down  toward  the  target, 

; rough  surface  I.  As  the  transducer  moves  across  the  surface  the  sound  pulse 
is  scattered  from  various  surface  features.  The  pictures  of  the  scattered 
field  (B,  C and  D)  change  as  the  surface  topography  changes.  In  each  case 
the  area  insonified  was  at  least  three  surface  correlation  lengths  in  extent. 
Insonifying  different  parts  of  the  surface  corresponds  to  sampling  a different 
population  (in  this  case  of  height  and  slope  distributions)  from  the  same 
statistical  ensemble.  The  predominant  features  in  all  visualizations  of  the 
scattered  field  are  the  angular  diffusion  of  the  sound  due  to  the  orientation 
of  the  surface's  slopes  and  the  interference  pattern  that  results  when  returns 
from  various  parts  of  the  surface  interfere  destructively  because  of  differing 
path  lengths.  The  scattered  field  in  Fig.  4 is  the  total  field  containing 
both  coherent  and  incoherent  components.  For  surface  I and  a 1 MHz  sound 
pulse,  the  coherent  term  which  includes  interference  dominates. 

Steady  State  Analysis. 

Experiments  in  microacoustics  have  traditionally  used  long  pulses  to 
make  measurements  at  the  resonant  frequency  of  the  source  transducer.  The 
signal  is  read  far  enough  into  the  pulse  (in  time)  to  make  transient  effects 
unimportant.  The  reflection  coefficient  for  the  rough  surface  is  determined 
by  normalizing  the  scattered  return  with  a return  from  a flat  surface  of  the 
same  material  at  the  same  distance  from  the  source  and  receiver.  Figure  5 shows 
reflection  coefficients  measured  at  two  frequencies  by  steady  state  techniques 
for  rough  surface  I.  The  geometry  shown  in  Fig.  2 is  retained  as  the  transducer 
is  moved  across  the  rough  surface.  Measurements  of  the  sound  scattered  back 
to  the  receiver  were  made  at  thirty-one  discrete  locations.  The  reflection 
coefficients  computed  for  these  points  are  joined  to  provide  the  solid  curves 
displayed  in  the  figure.  (At  normal  incidence,  scattering  back  at  the  source 
is  a case  of  specular  scattering  rather  than  of  backscattering. ) As  with 
schlieren  visualization,  there  is  no  means  to  separate  the  coherent  and  incoherent 
components  of  the  scattered  field.  However,  at  both  1 MHz  and  2.25  MHz  the 
coherent  component  dominates  the  scattered  field,  albeit  less  so  at  2.25  MHz. 

Transient  Pulse  Analysis. 

This  technique  makes  use  of  the  most  recent  advances  in  microacoustic 
I data  acquisition  and  processing.  A pulser  is  used  to  trigger  a transducer. 

This  generates  a short  (1.5-3. 5 cycle)  pulse  containing  a broad  spectrum  of 
frequencies  centered  on  the  resonant  frequency  of  the  transducer.  The 
experimental  procedure  used  consisted  of 

1.  sampling  the  waveform, 

2.  taking  a Fourier  transform,  and 

3.  computing  the  desired  reflection  coefficient. 

A series  of  31  measurements  were  taken  at  different  places  along  the  rough 
surface.  The  reflection  coefficients  were  computed  for  the  different  terms 
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(total  and  coherent)  In  the  following  ways: 

1.  Total  (includes  all  energy  regardless  of  path  length  or  phase) 

a.  Normalize  each  of  the  31  samples  by  the  proper  flat  surface 
return.  (Equivalent  to  computing  the  reflection  coefficient 
for  each  position  along  the  surface.) 

b.  Average  the  normalized  reflection  coefficients  and  determine 
the  mean  and  standard  deviation. 


[ Re  (Pgl)]  2 + [ Im  Cpsi> ] 2 

I Re  (Pref>l  2 + [I®  (Pref)l  2 


2.  Coherent  (retains  the  phase  term  through  the  averaging  process) 

a.  Average  the  real  and  imaginary  parts  of  the  scattered  returns 

separately.  (P  .) 

si 

b.  Compute  the  magnitude  of  this  average  scattered  return. 

c.  Normalize  the  result  by  the  magnitude  of  the  reference  pulse. 


where  Pgi  = Pressure  scattered  by  the  rough  surface  at  the  ith  position 

= Pressure  reflected  by  smooth  surface. 

A block  diagram  of  the  transient  pulse  acquisition  scheme  is  shown  in  Fig.  6. 

By  using  several  transducers  with  different  resonant  frequencies  (0.5,  1.0, 

2.25,  5.0  MHz)  it  was  possible  to  span  frequencies  from  0.14  MHz  to  8.0  MHz 
with  a resolution  of  0.025  MHz. 

Roughness  Parameter. 

Comparisons  in  rough  surface  scattering  are  made  through  the  use  of  a 
dimensionless  factor  known  as  the  Rayleigh  roughness  parameter  (referred  to 
herewith  as  roughness) . This  parameter  scales  the  roughness  according  to  the 
ratio  of  the  surface  rms  height  to  the  wavelength  of  the  sound.  This  is  then 
modified  by  the  angle  at  which  the  measurements  are  made.  The  form  of  the 
roughness  parameter  used  in  the  work  is 

2u 

Roughness  = — - a cos  0 = l 

A 1 

where  cos  0^  = 1 for  normal  incidence.  In  recent  years  measurements  have  been 
made  for  roughness  less  than  3.0;  however,  only  measurements  below  roughness 
of  1.5  were  in  agreement  with  theory.  Radar  measurements  characterized  the 
largely  incoherent  energy  for  very  large  roughness  (^  10.0).  Relative  influences 
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of  coherent  and  incoherent  terms  in  the  region  between  2.0  < E < 10.0  had  not 
been  explored.  Measurements  completed  through  this  task  covered  the  region 
from  E = 0.2  to  >_  20.0,  including  the  separation  of  the  coherent  component 
from  the  total  scattered  field. 

Results  and  Comparisons  with  Theory. 

Figure  7 shows  a summary  of  theory  and  experiment  in  terms  of  the 
reflection  coefficient  for  small  roughness.  Both  the  DeSanto  theory  and  the 
characteristic  function  involve  only  the  coherent  term.  This  is  also  true  of 
the  experimental  measurements  done  by  Clay,  Medwin  and  Wright.  The  solid  curve 
marked  as  the  total  reflection  coefficient  was  measured  in  this  experiment. 

Its  inclusion  in  this  illustration  indicates  a possible  difficulty  with  the 
Boyd  and  Deavenport  data.  Although  assumed  to  be  coherent,  the  data  falls 
closer  to  the  total  reflection  coefficient  curve  which  includes  both  coherent 
and  incoherent  terms. 

The  coherent  reflection  coefficients  measured  for  E < 3.0  are  shown  in 
Fig.  8.  The  smaller  values  of  roughness  («)  were  achieved  by  using  a 0.5  MHz 
transducer.  The  remaining  measurement  (x)  employed  a 5.0  MHz  transducer. 

These  measurements  continue  to  values  roughness  greater  than  20.  At  E = 1.9 
the  experimental  measurements  begin  to  differ  significantly  from  theoretically 
predicted  values.  Excursions  in  magnitude  that  begin  in  the  measurements  at 
E = 1.6  are  a characteristic  of  the  coherent  reflection  coefficient  as  the 
roughness  increases.  At  about  E = 2 the  interference  of  scattered  returns 
having  different  path  lengths  becomes  important.  The  measurements  shown  in 
Fig.  8 agree  with  theoretical  predictions  to  a far  greater  extent  than  the 
Boyd  and  Deavenport  data  in  Fig.  7.  Below  E = 1.5  the  measurements  should 
agree  with  theory.  Above  E = 2.0  where  interference  becomes  important,  neither 
theory  should  be  able  to  accurately  predict  the  behavior  of  the  sound  since 
neither  theory  takes  into  account  phase  differences  and  surface  correlation 
lengths  which  contribute  to  interference. 

The  measurements  shown  in  Fig.  8 are  extended  for  the  remainder  of  their 
range  in  Fig.  9.  The  presence  of  interference  minima  is  characteristic  for 
all  large  roughness,  coherent  measurements;  however,  the  location  and  extent 
of  the  minima  depend  upon  the  particular  surface  features  insonified.  In  other 
words,  the  scattered  field  for  different  statistical  samples  of  the  same 
ensemble  will  differ  in  some  detail  while  retaining  the  same  general  characteristics. 

The  total  reflection  coefficients  computed  from  the  same  measurements  of 
the  scattered  field  are  shown  in  Fig.  10.  The  average  value  of  the  reflection 
coefficient  is  bracketed  above  and  below  by  one  standard  deviation.  The  presence 
of  incoherent  terms  which  do  not  exhibit  phase  cancellation  virtually  eliminates 
any  noticeable  interference  pattern.  For  values  of  roughness  greater  than  2.0, 
the  incoherent  term  begins  to  dominate  the  reflection  until  at  E = 11.0,  the 
coherent  term  is  20  dB  down  from  the  total  reflection  coefficient. 

Figures  11  and  12  present  the  coherent  and  total  reflection  coefficients, 
respectively,  for  surface  I.  Since  the  rms  height  of  surface  I is  roughly  a 
factor  of  three  greater  than  for  surface  II,  the  range  of  roughness  measured 
with  the  same  acoustic  sources  is  a factor  of  three  larger  with  a corresponding 
decrease  in  resolution  (in  roughness,  not  in  frequency).  The  dominance  of  the 
incoherent  term  at  large  roughness  is  evident  in  this  case  also.  Direct 
comparisons  can  be  made  between  the  total  reflection  coefficients  for  both 
surfaces  using  the  roughness  parameter  as  the  appropriate  scale.  The  same  is 
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not  true  for  the  coherent  reflection  coefficients  at  S >1.5.  This  is  the 
region  in  which  path  length  interference  begins  to  occur.  Because  individual 
surface  features  become  important,  the  surface  correlation  length  enters  as 
a necessary  factor  in  the  comparison.  If  the  roughness  parameter  is  further 
normalized  by  surface  correlation  length,  the  coherent  reflection  coefficients 
agree  even  for  large  roughness. 

Summary. 

The  NRL  experiment  has  (a)  begun  in  the  low  frequency  region  where 
comparisons  with  the  characteristic  curve  would  provide  a validation  of  the 
measurement  scheme,  (b)  spanned  the  region  of  intermediate  roughness  where 
previous  data  were  questionable,  and  (c)  provided  scattering  measurements  at 
high  roughness  where  previous  data  were  inadequate.  Both  the  total  reflection 
coefficient  and  the  coherent  term  were  computed  from  the  same  set  of  measure- 
ments. Coherent  coefficients  are  needed  for  comparison  with  many  of  the  newer 
theoretical  models  and  for  assessing  the  relative  importance  of  height  and 
slope  distributions,  correlation  lengths  and  source  characteristics  as  controlling 
parameters  in  the  scattering  process.  The  total  reflection  coefficient  is  needed 
for  application  to  a variety  of  practical  problems,  particularly  those  concerned 
with  sound  propagation  in  the  Arctic  environment. 

Work  Distribution  Between  NRL  and  A.U. 

All  experimental  work  was  done  at  NRL  together  with  the  data  analysis 
and  reporting.  Computer  software  and  hardware  developed  through  the  contract 
at  The  American  University  made  possible  the  comparisons  with  theoretical  models 
and  implemented  the  experimental  design  used  in  the  transient  pulse  analysis 
measurements. 
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TABLE  1.  PARAMETERS  OF  ROUGH  SURFACES  I AND  II 


Parameters 

Surface  I 

Surface  II 

RMS  height  (0) 

0.124  CM 

0.037  CM 

Correlation  Length 

0.737  CM 

0.305  CM 

Maximum  Gradiant 

0.492 

0.458 

Maximum  Height 

0.238  CM 

0.094  CM 

Maximum  Depth 

-0.313  CM 

-0.078  CM 

Minimum  Radius  of  Curvature 

0.891  CM 

0.328  CM 

Milling  Accuracy 

0.005  CM 

0.005  CM 

Mean-Squared  Deviation  from  the 
Theoretical  Distribution 

0.0036 

0.0016 

Photograph  of  Surface  I.  The  rms  height  is  0.124  cm  with  a 0.737  cm  correlation 
length.  The  block  is  25.4  x 7.6  x 7.6  cm.  THe  milled  rough  surface  is  17.0  cm  long 
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GEOMETRY  OF  THE  EXPERIMENT 


Fig.  2 The  orientation  of  a source  and  surface  for  all  the  experi- 
ments is  pictured.  The  sound  was  normally  incident  on  the 
nominal  mean  of  the  surface. 
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Fig.  3 A.  typical  schlieren  visualization  system  is  shown  as  a block 

diagram.  Disturbances  in  the  water  caused  by  the  sound  pressure 
produce  an  image  at  the  camera  lens. 


Schlieren  photographs  of  the  incident  and  the  scattered  fields.  The  incident  pulse 
is  shown  in  (a).  The  scattered  field  in  (b),  (c),and  (d)  demonstrates  the  effects 
of  phase  cancellation. 


TRACING  OF  RANDOM  SURFACE 
CORRELATION  LENGTH' '0.74  CM 


The  reflection  coefficients  are  shown  above  the 


transducer  diameters  and  the  surface  topography  (surface  I) 


are  drawn  to  scale 
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RAYLEIGH  ROUGHNESS  PARAMETER 

Fig.  8 The  measured  coherent  reflection  coefficients  for  surface  II 
are  plotted  against  the  Rayleigh  roughness  parameter. 
Measurements  made  with  the  0.5  MHz  transducer  are  indicated 
by  solid  dots  (•) . The  5.0  MHz  transr'-jcer  was  used  for  the 
remaining  data  (X) . The  characteristic  curve  (dashed  curve 

) and  DeSanto's  model  (solid  curve  — ) are  shown  for 

comparison  as  in  Fig.  7. 
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Fig.  10  Total  reflection  coefficients  for  5.0  MHz  transducer  obtained 


by  transient  pulse  analysis  are  plotted  against  the  Rayleigh 


roughness  parameter  for  brass  surface  II.  The  standard  deviation 


for  each  data  point  brackets  the  mean  curve 
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